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Abstract 

Starting from the second law of thermodynamics applied to an isolated system consisting of 
the system surrounded by an extremely large medium, we formulate a general non-equilibrium 
thermodynamic description of the system when it is out of equilibrium. We then apply it to 
study the structural relaxation in glasses and establish the phenomenology behind the concept of 
the fictive temperature and of the empirical Tool-Narayanaswamy equation on firmer theoretical 
foundation. 
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I. INTRODUCTION 



It is well known that when a liquid is disturbed suddenly from its equilibrium state by 
changing the temperature or pressure of the surrounding medium, or both, then the liquid 
undergoes a rapid, solidlike change, followed by a slower, liquidlike change towards the 
new equilibrium state. These changes can be seen in the variation in its thermodynamic 
properties such as the volume V or its enthalpy H with time. For a supercooled liquid, the 
above scenario plays an important role. As the temperature is lowered or the pressure is 
increased, the scale separation between the fast and the slow processes in supercooled liquids 
increases until the latter becomes too large compared to the experimental observation time 
r obs . In this case, the system is said to be kinetically arrested in that the liquidlike changes 
no longer contribute to the observed properties. The system behaves like a solid and is 
called a glass Jfl S, fl Q. 

The glass is a system that may be far from equilibrium so one cannot apply equilibrium 
statistical mechanics to investigate its properties, which vary with time. One must resort to 

nn 

apply non-equilibrium thermodynamics [3,[8(, not a well-established field at present, to study 

ly known as structural relaxation. It 
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glasses and their relaxation in time; the latter are usua 
is known from the early work of Littleton jt| , and Lillie 
glass changes during relaxation, thus implying the dependence of the relaxation on changes 
in the state of the glass The most general framework for developing non-equilibrium 
thermodynamics must satisfy the second law of thermodynamics or must start from it. The 
generality arises since the second law is independent of the processes (structural or otherwise) 
going inside the system. It is also independent of the details of the system considered and 
does not requires any sophisticated concepts like ergodicity or its loss, etc. Our main goal in 
this work is to develop an appropriate non-equilibrium thermodynamics, which will then be 
applied to glasses with the hope to gain some new insight and to clarify at the fundamental 
level certain concepts extensively used in glasses. 

Glassy behavior and their properties have been extensively studied and usually explained 
by invoking empirical rules 6, Q] that, although they have proved invaluable and very 



reliable, lack theoretical justifications |3|,|4]. Only recently, attempts have been made [l4, Il5] 
that use modern framework to investigate non-equilibrium properties of glassy relaxation at 
the molecular level. Our treatment here differs from these attempts in that we develop our 
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approach using the second law that should be applicable to all systems including glasses. We 
do not derive the actual laws of relaxation for which one must turn to other sources such as 
[lil [l5|. Our goal is quite different. We wish to understand some of the important concepts 
used for glasses at a fundamental level. One of the most widely used concept in this field 
is that of the ftctive temperature, first introduced by Tool [3] in an empirical fashion to 
describe non-linear relaxation in glasses. The system under study slows down so much upon 
reducing the temperature from its initial value T', where the system was in equilibrium, 
to some temperature T that one has to wait for a very long time before true equilibrium 
is reached at the final temperature. In this case, crudely speaking, the glass properties 
are assumed to be similar to a fictive liquid at some intermediate temperature between T' 
and T . As time goes on and as the system undergoes structural rearrangements to come to 
equilibrium, the temperature of the system continues to change and finally becomes To. This 
already means that the fictive temperature of the system continuously changes from T' to To. 
Despite its continual usage in the field, the true meaning of the fictive temperature, though 
phenomenologically obvious, is not well-defined in terms of fundamental quantities such as 
the entropy. In particular, there exists a variety of fictive temperatures, each associated 
with the relaxing quantity under investigation, which makes the concept not very rich. 

There is another aspect of structural relaxation. Its presence means that the glass is a 
non-equilibrium state. Thus, its temperature must be changing during the process of relax- 
ation. How does one define the instantaneous temperature of the liquid? The instantaneous 
temperature itself must relax to To as time goes on. Thus, there will a relaxation time de- 
scribing the relaxation of the temperature of the glass. Tool [12] and Narayanaswamy 13] , 
among others, observed that the relaxation time not only depends upon the temperature 

T), but also depends upon the fictive temperature of the system; see for example qs 3. is 
the instantaneous temperature of the glass the same as the fictive temperature? These are 
important issues as a deeper understanding of these concepts will provide a more qualitative 
and predictive understanding of glass transition. 

The layout of the paper is as follows. We consider an isolated system consisting of the 
system of interest surrounded by a very large medium and follow the consequences of it in 
the next section. In Sect. 1111} we follow the consequence of partial equilibrium to develop 
a very general non-equilibrium thermodynamics, which is then applied to a glass in Sect. 
IIVI The concept of the fictive temperature and the Tool-Narayanaswamy phenomenology 



3 



are considered in Sect. |Vj and established on a firm theoretical ground. The conclusions are 
given in the last section. 



II. CONSEQUENCES OF THE SECOND LAW 

As said above, we study non-equilibrium systems by proceeding in a general manner by 
following the consequences of the second law, which is well established. As usual, we apply 
the second law to an isolated system, which we denote by £ ; it consists of the system E 
of interest (such as our glass) in a medium denoted by S containing it. We will consider a 
single component system, which is sufficient for our purpose. According to the second law, 
the entropy So of an isolated system S can never decrease in time [if] ]: 

dSn(i) , . 

-^ 1 >0. 1 
dt ~ w 

What happens inside the isolated system (loss of ergodicity in parts of the system, chemical 
reactions, phase changes, etc.) cannot affect the direction of the inequality, which makes 
it the most general principle of non-equilibrium thermodynamics. The law itself imposes 
no restriction on the actual rate of entropy change. In general, So also depends on the 
number of particles N , energy E , and volume V of S . Thus, S (t) used above should be 
really written as Sq(Eq, Vq, No, t). However, as the extensive quantities remain constant in 
time there is no harm in using the compact form So(t) during approach to equilibrium. The 
entropy So(t) is a continuous function of each of its arguments. The energy, volume and 
the number of particles of S are denoted by E, V, and N, respectively, while that of the 
medium E by E, V, and N. Obviously, 

E = E + E, V = V + V, N = N + N. 

We will assume that N of the system is also fixed, which means that N is also fixed. However, 
the energy and volume of the system may change with t. 

When the isolated system is in equilibrium, its entropy Sq(E , Vq, N , i) has reached its 
maximum and no longer has any explicit time-dependence so that it can be simply written 
as So(E , Vo, Nq) or So- In this case, different parts of S have the same temperature T and 
pressure P : 

1 _ dSo Pq _ dSo , s 

To ~ dEo' To ~ Wo' U 
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we have defined the temperature by setting the Boltzmann constant ks = 1 in this work. 
Otherwise, the entropy So(t) continuously increases and the isolated system is said to be 
not in equilibrium. The medium is considered to be very large compared to £, so that its 
temperature, pressure, etc. are not affected by the system. We assume £ to be in internal 
equilibrium (its different parts have the same temperature and pressure, but £ and £ may 
not be in equilibrium with each other). Thus, its entropy S no longer has an explicit time 
dependence, but has an implicit t-dependence through the t-dependence of E, and V. The 
time variation of S (t) is due to the relaxation going on inside £ as it is driven towards 
equilibrium with the medium. 

The entropy S (t) of the isolated system can be written as the sum of the entropies S(t) 
of the system and S(t) of the medium: 



S (E , Vo, N , t) = S(E, V, N, t) + S(E, V, N); 



(3) 



there is no explicit t-dependence in S(E, V, N) due to internal equilibrium. The correction 
to this entropy due to the weak stochastic interactions between the system and the medium 
has been neglected, which is a common practice [lfjj. We expand So in terms of the small 
quantities of the system 1(| 



S(E,V,N)~S(E ,V ,N)- ( || 



i dS 

m - 1 w 



V(t). 



Vo 



It follows from the internal equilibrium of £ that 



OS 
dE 



to 



1 



dS_ 
dV 



Vo 



Po 
To' 



and S = S(Eq, V ,N), which is a constant, is independent of the system. Thus, 



S (t) -S~ S(E, V, N, t) - E(t)/T - P V(t)/T . 



(4) 



Let us introduce 



G(t) = H(t) - T S(t), H(t) = E(t) + PoV{t), 



(5) 



the time-dependent Gibbs free energy and enthalpy of the system £ with the medium £ at 
fixed Tq and Pq. We thus finally have 



S (t) -S = S(t) - H{t)/To = -G(t)/T 



(6) 



so that the behavior (pQ) of So(t) of the isolated system leads to an very important conclusion 
about the Gibbs free energy of the system: 

dG(t) 

The Gibbs free energy G(t) decreases as the system relaxes towards equilibrium, a result 
quite well known. If we abruptly cool the system from some previous temperature such as T g 
to a lower temperature, the Gibbs free energy at the new temperature remains equal to its 
value at the previous temperature at time t — 0. As it relaxes, G(t) continuously decreases. 
It cannot increase without violating the second law. 

We have given the essential steps in its derivation here not only for the sake of continuity 
as some of the intermediate steps will be needed later on, but also to make some important 
points, which we now list. 

1. In deriving the above equation ([6]), no assumption about the system E has been made. 
In particular, we have not assumed any particular aspect of its non-equilibrium nature, 
such as a particular form of relaxation (Arrhenius or otherwise), loss of ergodicity, etc. 

2. The identification of So(t) — S with the Gibbs free energy G(t) of £ is generally valid 
under the assumption of the medium being large compared to S, which can be satisfied 
as well as we wish. 

3. The Gibbs free energy G(t) and the enthalpy H(t) are determined by the temperature 
T and the pressure Po of the large medium. 

4. The continuity of So(t) with respect to all of its arguments that was mentioned earlier 
also applies to the Gibbs free energy of the system. 

5. The decrease in G(t) must not be violated even when there is a loss of ergodicity in 
the system, as is commonly believed to occur in a glass transition. 

6. For glasses, we have an additional experimental fact. The enthalpy remains continuous 
across the glass transition. The continuity of the enthalpy then implies that the entropy 
S(t) also will remain continuous with respect to all of its arguments, as the system 
relaxes. This allows us to differentiate the entropy, which will be required in Sect. IIHI 
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III. NON-EQUILIBRIUM THERMODYNAMICS 



When the equality in ([I]) occurs, different parts of So (such as E and E) have the same 
temperature T and pressure P . Otherwise, they have different temperatures and pressures, 
in which case a common assumption made b y a lmost all workers is that of partial equilib- 
rium (see, for example, Landau and Lifshitz ly, see p. 13]) when Eo is out of equilibrium; 
each part is in internal equilibrium (local equilibrium), which then allows us to define the 
temperature, pressure, etc. for each part, which may all be different. In this situation, their 
entropies have the maximum possible values for their respective energies and volumes, and 
the number of particles. As a result, they have no explicit t-dependence [see the equilibrium 
condition (|2J) above for E ]; their variation in time comes from the time variation of their en- 
ergies, volumes, etc. The entropy S of the system determines its instantaneous temperature 
T{t) and pressure P(t): 

ds_ j_ oi = m 

dE T(tY dV T(t)' { ' 

These are standard relations for the entropy [161 ]. except that all quantities except S in 
the above equations may have an explicit dependence on time t that will make S depend 
implicitly on time. Accordingly, 

dS 

under internal equilibrium. Relations like (jHJ) along with (J9]) for internal equilibrium are 
used commonly in non-equilibrium thermodynamics. For example, we use them to establish 
that heat flows from a hot body to a cold body; see for example Sect. 9 in Landau and 
Lifshitz [16J. The glassy state, in which the fast dynamics has equilibrated and the slow 
dynamics is extremely slow, will thus be treated as a state in internal equilibrium, although 
it is not in equilibrium (with the medium). This observation will be very important when 
we discuss the concept of the Active temperature in Sect. [Vj 

Recognizing that S(t) has no explicit t-dependence, see (jSj), but is a function of E{t) and 
V(t) (N is kept a constant), we have for the differential dS(t) 

dS(t) = ^- ) dE(t) + ^ ) dV(t), 

where we have used (JS} and have allowed the pressure and the temperature of the sys- 
tem to be different from those of the medium for the sake of generality. The first law of 
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thermodynamics follows from this equation: 

dE(t) = T(t)dS(t) - P(t)dV(t), (10) 

which does not depend on the temperature and pressure of the medium. In this form, the 
first law has the standard look with the first term representing the heat 

dQ = T(t)dS{t) (11) 

added to the system and the second term without the sign denoting the work 

dW = P(t)dV(t) 

done by the system. 

Using H(t) = E(t) + P V(t) in ®, we find that 

dH(t) = T(t)dS(t) + V(t)dP + [P - P(t)]dV(t), (12) 

where the last term appears due to the lack of equilibrium with the medium. Accordingly, 
the heat dQ, see (ITTil . is no longer equal to dH{t) at constant pressure Pq of the medium: 

dQ(t) = dH(t)\ Po + [P(t) - P ]dV(t). 

The specific heat Cp at constant pressure is given by 

However, obtaining the entropy of the system from the measured values of the specific heat 
requires care: 

dS{t) = Tj^dTo < jr dT Qi 
which follows from ( Tl9l) . derived later. 

The differential of G{t), see (jSJ), turns out to be 

dG(t) = -S(t)dT + V{t)dP + [T(t) - T )dS(t) + [P - P(t))dV(t), (13) 

Again, the last two terms are corrections to dG(t) due to the non-equilibrium nature of the 
system. We observe from (ITS]) that 
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Again, the last two terms in each equation are the correction due to non-equilibrium nature 
of the process, and would be absent in an equilibrium process. 

One can compare the Gibbs free energy differential in (TT5|) with the approach developed 
by de Donder 17] and Prigogine The last two terms in (1131) denote the contributions 
from two different structural order parameter or the degree of advancement. In the present 
context, the two parameters are determined by the instantaneous entropy and volume 

_ S(t)-S(oo) _ V(t) - V(oo) 
4l " S{0) - S(oo) ' 42 " V{0) - V(oo) ' 

and the corresponding affinities are given by 

Al = _[ T ( t ) _ To ] [5(0 ) _ 5(00)] , A 2 = [P(t) - P ] [V(0) - V(oo)] . 
Thus, we have 

[T(t) - T ]dS(t) - [P(t) - P ]dV(t) = - A 2 £ 2 , 

so that each contribution A^ < 0, as expected from the variation of G(t) during relaxation 
at constant To, Pq. We also note that we can write the first law as 

dE{t) = T dS(t) - P dV(t) - A& - A 2 &, (14) 

as expected from the standard formulation by de Donder. 



IV. RELAXATION BELOW THE GLASS TRANSITION 

We now apply the above formalism to a glass. Above the glass transition temperature 
T g but below the melting temperature, the system is a supercooled liquid (SCL) as the 
relaxation time r of S remains less than the observation time r b s . At T g , they become 
identical. Below T g , r becomes larger than r b s , and the system turns into a glass. Let us 
consider the system in the glassy state. With time, the glass (E) will relax so as to come 
to equilibrium (the corresponding SCL state obtained by increasing r b s to the relaxation 
time at that temperature) with the medium if we wait longer than r b s . It should be noted 
again that, due to the internal equilibrium of the system, there is no explicit t-dependence 
in S on the right side of (EJ) or (@J. Accordingly, the t-dependence in S, H , and G is implicit 
through E(t), and V(t). In the following, the glass is considered to be formed under isobaric 
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conditions so the pressure of the medium is kept fixed at Po as its temperature is step-wise 
varied. Accordingly, its instantaneous pressure P(t) is always equal to Po of the medium 

Pit) = P , (15) 

but its temperature will in general be different than To and vary in time, as we will show 
below. The initial enthalpy H(0) at T is the enthalpy of the glass at temperature T', and 
H(oo) the value of the SCL enthalpy after complete relaxation at temperature T . It is 
experimentally found that the enthalpy decreases with time during an isothermal relaxation 
so that 

H(t = 0) > H(t -> oo). (16) 

This decrease is a general property of thermodynamics which follows from the specific heat 
being non-negative. We first consider the case when T' is the glass transition temperature, 
so that the system is in the SCL state at T' . As T' > T , the enthalpy H(0) of SCL at T' 
must be higher than the enthalpy Hit — > oo) of SCL at T . The same is also true of the 
volume in many cases, which relaxes to a smaller value in an isothermal relaxation. However, 
this property of the volume is not a thermodynamic requirement. Accordingly, as a general 
rule 

dH(t) , . 

during isothermal structural relaxation in glasses. In the following, we will only use the 
above general property of the enthalpy, and not of the volume. Now, if T 1 is below the glass 
transition, then from the result just derived, we conclude that H(0) is even larger than the 
SCL enthalpy at T' . This even strengthens the above inequality (fT6j) . 



Let us now turn to the time derivative of the entropy Sq, which is changing because the 

energy and volume of S are changing with time [if]]. Thus, 

dS Q (t) _dS__ 1 dEjt) _ P dVjt) 
dt dt To dt Tq dt 

" \dE T J dt \dV T J dt ~ ' 
as the relaxation goes on in the system E. It is clear that 

dS_ , 1 OS P 

if dSo/dt > 0. Thus, as long as the relaxation is going on due to the absence of equilibrium, 
the two inequalities must hold true. Accordingly, the derivative dS/dE, which by definition 
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FIG. 1: Schematic behavior of the entropy for SCL (blue curve) and GL (red dotted curve). The 
GL entropy decreases, shown by the downward arrow, as it isothermally (constant temperature Tq 
of the medium) relaxes towards SCL, during which its temperature T(t) also decreases towards Tq. 



represents the inverse temperature 1/T{f) of the system, see (ED, must be different from 
1/Tq of the medium: 

T{t) ± T . 

As dS/dV = Po/T(t), we see immediately that 

dS (t) ( 1 1\ dH(t) 

From ffTTl) . we observe dH(t)/dt < during relaxation in glasses. Thus, we are forced to 
conclude that 

T(t) > T , (19) 

the equality occurring only when equilibrium has been achieved. The instantaneous tem- 
perature can be measured by using a small "thermometer" so as not to disturb the internal 
equilibrium of the glass. Such a measurement will allow us to explore its variation in time. 
The above calculation also shows that 

dS(t) 1 dH(t) 
dt T(t) dt ' 1 J 

which is the first term in (|T8|) . The equation above can also be obtained from (|T2l) . Using 

P(t) = Pq for an isobaric process, which is the normal situation in most experiments, we see 
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that the last term in (fT2l vanishes. Hence, 



dQ = dH{t) = T(t)dS(t) 

is valid in all isobaric processes, from which (|20|) follows immediately. 

The relaxation that occurs in the glass originates from its tendency to come to thermal 
equilibrium during which its temperature T(t) varies with time; recall that we are considering 
a cooling experiment. The relaxation process results in the lowering of the corresponding 
Gibbs free energy, as is seen from (jjj), which is a consequence of the second law in (CD). 
Accordingly, there are changes in its enthalpy and entropy, which are in the same direction; 
see (j20p . The lowering of G(t) with time results in not only lowering the enthalpy in a 
cooling experiment, as observed experimentally, but also the entropy S(t) during relaxation: 

(dS(t)/dt) < 0, 

as shown in Fig. [IJ At constant P an d T , we see from (fT3|) that 

dG(t) = [T(t)-T ]dS{t), (21) 

from which it follows that 

showing that G(t) converges to its equilibrium value more slowly compared to the conver- 
gence of S(t) as T(t) — > To, i.e. as t — > oo as the above derivative vanishes in this limit. 

It is instructive to compare the specific heat of the glass with the specific heat of the 
corresponding fully relaxed state obtained as t — > oo. Let us assume that at time t — 0, we 
change the temperature of £ form some initial temperature T' to To instantaneously. We 
consider the system at t — t Q b s and determine its enthalpy. The specific heat of the glassy 
sample at this instant is given by 

Cpg= , im g<°)-gfa.). 

' g AT^O AT 

Then, the corresponding specific heat after complete relaxation is given by 

Cp relax = Hm " — ~ . r-r—i — " ~ ^ Cp e . 

' AT^O AT ~ ,g 
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V. THE FICTIVE TEMPERATURE AND THE TOOL-NARAYANASWAMY 
EQUATION 



We now need to turn our attention to the distinction between the fast and slow degrees 
of freedom (dof), a characteristic of any glass. Situation similar to this also occurs in the 
attainment of thermal equilibrium between the nuclear spins and their environment during 
nuclear relaxation 19( , where the spin-lattice relaxation is extremely slow. The explanation 
of this kind of behavior (slow and fast dof) in a wide class of substances lies in internal 
molecular motions other than simple vibrations. The fast dof cool down and equilibrate 
very fast, while the slow dof take much longer to transfer their energy and equilibrate 
because of very weak coupling with the surrounding medium. Here, we are talking about 
equilibration with the medium. We will denote those dof that have equilibrated with the 
medium at time t by a subscript "e", and the remaining that are not equilibrated by "n." 

We set t = at the instant the system is abruptly cooled from its equilibrium state at T' 
to some lower temperature T of the medium. At this time, all dof are out of equilibrium 
with the medium at To. The fast dof equilibrate within the observation time t b s , with the 
slow dof remaining out of equilibrium [20]. Eventually, all dof come to equilibrium with the 
medium. Thus, the number of dof in equilibrium keeps on increasing with time. Let D 
denote the total number of the dof in the system, which is determined by the number of 
particles N in it; hence, it remains constant. Let D e (t) and D n (t) denote its partition in 
equilibrated and non-equilibrated dof, respectively: 

D = D c (t) + D n (t); 

evidently, they keep varying in time. As said above, the clear distinction between the two 
kinds of dof arises because of a very weak coupling between them and of the slow dof with 
the medium. The weak coupling allows us to treat them as almost uncorrelated and quasi- 
independent, which then immediately leads to the following partition of the entropy, the 
energy and the volume into two contributions, one from each kind because of their quasi- 
independence mentioned above: 

S(t) = S e (t)+S a (t), (22a) 
E{t) = E e {t) + E n {t), (22b) 
V(t) = V c (t) + V n (t), (22c) 
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where the notation is self-evident. It should be noted that S e (t) and S n (t) stand for 
S c (E e (t), V e (t)) and S n (E n (t), V n (t)). The corrections to each of the partition due to the 
very weak coupling is small enough to be neglected. 

One should not confuse dof c with only fast dof. To see this most clearly, we recall that at 
t = 0, none of the dof have equilibrated at the new temperature T , so that D c (t = 0) = 0. 
But all of the same dof were equilibrated at T", implying that dof e at T' contains fast and 
slow dof. The same happens as t — > oo, in which case D e (t) — > D at To implying that all 
dof, fast and slow, have equilibrated. Thus, in general, D e (t) contains both fast and slow 
dof. Let us now consider dof n . At t = 0, none of the dof have equilibrated at the new 
temperature T . Accordingly, D n (t = 0) = D, so that dof n contains both fast and slow dof. 
However, for t > t Q b s , only (or mostly) the slow dof remain in D n (t). 

Let us now introduce the following derivatives of the energy partitions 



dEJt) , , dEJt) , , 

^■m^-^sm- (23) 



at a given t, so that 



9Se{t) - x(t) dSe{t) 9S ^-\1 x (t)] dS » {t) (24) 

dE(t) ~ [ ) dE e (ty dE(t)~ [ [ n dE Q (ty [ ] 

The derivatives in the two equations above are at fixed V e (t) and V n (t), respectively. At 
t = 0, D e {t = 0) = 0, E e {t = 0) = and x(t = 0) = 0. At t -»• oo, D c (t) -> D, E n (t) = so 
that x(t) = 1. As time goes on, more and more of the "n" dof equilibrate, thus increasing 
D e (t) and x(t). 

By definition, we have 

dS e (t) J_ 
dE e (t) T ' 

which follows from the equilibrium of the dof e with the medium, while the dof n will have a 
temperature different from this. Accordingly, we introduce a new temperature T n (t), defined 
by the derivative 

dE n (t) T n (ty { ) 

The following identity 

1 _x(t) + l-x® 



T(t) T T n (t) 

easily follows from considering dS(t)/dE(t) and using (I22al) and (|2^|) . Initially, x(0) = 
so that T(0) = T n (0) = T', while T(t) — >• T as t — >• oo, as expected. This division of 
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the instantaneous temperature T(t) into T and T n (t) is identical in form to that suggested 
by Narayanaswamy [13|, except that we have given thermodynamic definitions of the non- 
linearity parameter x(t) in (I23p and of the new temperature T n (t) (125]) in our approach. 
Both these quantities, being intensive, can only depend on energy and volume per particle, 
through which these quantities gain their implicit t-dependence. 

Let us now understand the significance of the above analysis. The partition of the ther- 
modynamic quantities in (1221) along with the definition of the fraction x(t) shows that the 
partition satisfies a lever rule: the relaxing glass can be conceptually (but not physically) 
thought of as a "mixture" consisting of two different "components" corresponding to dof e 
and dof n : the former is at temperature T and has a weight x(t); the latter with a comple- 
mentary weight 1 — x(t) is at a temperature T n (t). Thinking of a system conceptually as 
a "mixture" of two "components" is quite common inn theoretical physics. One common 
example is that of a superfluid, which can be thought of as a "mixture" of a normal vis- 
cous "component" and a superfluid "component" 2l|, Sect. 23]. In reality, there exist two 
simultaneous motions 21], one of which is "normal" and the other one is "superfluid". A 
similar division can also be carried out in a superconductor: the total current is a sum of a 
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Sect. 44]. 



"normal: current and a "superconducting current" 

The division of the dof envision above is no different from these divisions in a superfluid or 
a superconductor. However, because of the non-equilibrium nature of the system, there is an 
important difference here compared to a superfluid or a superconductor. The e-component 
is in equilibrium (with the medium), but the n-component is only in internal equilibrium. 
While the significance of the former as a SCL "component" at Tq,Po (dof=_D e ) is obvious, 
the significance of the latter requires clarification. At £ = 0, T n (t) represents the temperature 
T' of the equilibrated SCL (dof=D) from which the current glass is obtained by cooling. At 
this time, the entropy S n (t = 0) of initial state of the glass at T is equal to the entropy 
Sscl(T') of the equilibrated SCL (dof=D) at the previous temperature T'. The latter has 
the energy and volume Escl(T') — E n (t — 0) and volume Vscl(2~") = V n (t = 0). At any later 
time t > 0, T n (t) represents the temperature associated with the energy E n (t) and volume 
V Q (t) of the non-equilibrated "component" of the glass and has a weight 1 — x(t). This 
"component," being in internal equilibrium, can be identified as a fictive SCL [dof=Z) n (t)] 
at temperature T(t) < T' of energy -Escl = E Q (t) and volume Vscl = V^t). In other words, 
the relaxing glass at any time t can be considered as consisting of two SCL "components," 
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one at temperature T [dof=D e (t)} and the other one [dof=D n (t)} at temperature T n (t). 
The temperature T = T n (t) uniquely determines the energy Escl(T) = E Q (t) and volume 
^scl(^) = V n (t) of the corresponding fictive SCL [doi=D Q (t)}. 

As the fictive liquid at T = T n (t) contains only (or mostly) the slow dof, it does not 
yet really represent a SCL associated with the system at T = T n (t), as the former lacks 
dof e , while the latter contains all dof. This does not pose any problem as the missing 
dof e at T = T n (t) are in equilibrium not only with the dof n at T = T n (t), the fictive SCL 
mentioned above, but also with the medium at T = T n (t). Thus, one can consider "adding" 
these missing dof e (dof=_D e ) to the fictive liquid, which now represents the equilibrated SCL 
(dof=-D) at T = T n (t). This SCL is not the same as the glass with its fictive T n (t), as 
the latter has its dof e at T while the SCL has all of its dof at T Q (t). However, all of their 
thermodynamic properties associated with dof n must be the same, as their entropy function 
is the same for both liquids. Similarly, the SCL "component" at T ,P (dof=D e ) should 
also be "supplemented" by the missing dof n to give rise to the equilibrated SCL at T ,P 
(dof=£>). 

We are now in a position to decide which of the temperatures T(t) and T n (t) qualifies as 
the fictive temperature. This temperature is supposed to characterize the non-equilibrium 
aspect of the system. As T(t) contains information about both kinds of dof, it is not the 
appropriate temperature to be identified as the fictive temperature, even though it depends 
on t. The temperature T n (t), on the other hand, depends only on non-equilibrated dof n , 
and should be identified as the fictive temperature of the relaxing glass at time t. This 
temperature is not the internal temperature of the glass at this time, but represents the 
equilibrium temperature of the corresponding SCL at T = T n (t) as noted above. 

As first pointed out by Littleton |9( and Lillie lOj, and discussed by several authors, see 
for example 0, 0, 11], the viscosity keeps changing with time during relaxation. Thus, if one 



uses an Arrhenius form for the viscosity, it must depend not on To, but on T(t); it is the 
instantaneous temperature that characterizes the instantaneous state of the glass. Thus, the 
Arrhenius form for the viscosity, which is usually taken to be proportional to the relaxation 
time, must be expressed as 

B 



r)(t) = r] exp 



[T(t) 



7? exp 



x(t) 1 - x(t) 



(27) 



the form conventionally identified as the phenomenological Tool-Narayanaswamy form. 
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Here, 770 and B are some parameters of the system and may depend on To, Pq and also 
weakly on time t. Our derivation above justifies this form on a solid theoretical ground. 

One can carry out a similar analysis with decomposing the volume; see (j22cj) . However, 
we do not obtain any new result as Pit) — Pq- To see this, we proceed exactly as above but 
use the volumes instead of the energies. Introducing the parameter x v (t) defined by 

v[ ) - dv{ty 

at fixed E e (t), which may be different from x(t), and using 

dS e (t) P dS n (t) _ P 



dV e (t) T ' dV n (t) T^ v (ty 
at fixed E c (t) and E n (t), respectively, along with dS(t)/dV(t) = P /T(t), see © and (fT5l). 
we find the following decomposition of the inverse instantaneous temperature 

1 _ x v (t) 1 - x v (t) 



T{t) T T n>v (t) 
Now, the new Active temperature T n>v (t) represents the temperature T = T n ^(t) of corre- 
sponding Active SCL [dof= D n (t)] with energy and volume E SCL (T) = E n {t) and volume 
Vscl(T) = V n (t). This Active liquid is the same as noted above as far as the energy and 
volume are concerned. However, as SCL is an equilibrated state, the specification of energy 
and volume uniquely determines the temperature, which from the earlier analysis was seen 
to be exactly T n (t). Thus, we conclude 

T n , v {t) EE T n {t). 

In order for the above decomposition to be consistent with fT2"6T) at all times, we must ensure 
that 

x v {t) = xit). 

In other words, our definition of the Active temperature gives the same value whether we 
consider the energy relaxation or the volume relaxation. 

It is highly likely that the slow relaxation consists of many difierent relaxation modes, 
which we index by j = 1, 2, • • ■ . However, there does not seem to be any strong argument 
to suggest that all these difierent relaxation modes are almost decoupled as was the case for 
the fast and slow relaxations studied above. In that case, it is not possible to partition 
the thermodynamic quantities such as the entropy, etc. associated with dof n as a sum over 
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these different modes. Despite this, let us follow the consequences of such an assumption. 
We express D n {t) as a sum over j 

D n (t) = Y;D<p{t), 

j 

where the notation is quite transparent. We similarly express all of the n-quantities in ([22]) 
as a sum over j. We can similarly express 

1 -x(t) = ^Vj, 
j 

where 

_ dE$\t) 
Vj ~ dE(t) ' 

We can now introduce a Active temperature for each j-th n-dof 

dsP(t) 1 
dE%\t) T^\t) 

such that 

1 - x(t) = \- yj 

a decomposition also described by Narayanaswamy [13(. However, because of the above 
equality, the presence of more than one kind of relaxation modes does not change the earlier 
decomposition (1261) . In other words, no new insight is gained by such an assumption. One 
can introduce an equilibrated Active SCL at each of the Active temperatures as above. We 
will not stop here to do so as it is straightforward. 

VI. CONCLUSIONS 

We have developed a non-equilibrium thermodynamics to study systems away from equi- 
librium. The approach is quite general and is not limited to systems close to equilibrium. 
Assuming internal equilibrium, a common practice in the field, we find the correction to 
the differential free energies that are consistent with de Donder-Prigogine approach to non- 
equilibrium thermodynamics; see for example (114p . Even though, we have mainly discussed 
supercooled liquids, the approach does not require the presence of a melting transition and 
an equilibrium crystal for its application. Thus, it should also be applicable to other glassy 
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systems such as spin glasses, where there is no equilibrium crystal as the true equilibrium 
state. The only requirement is that enough time has passed after the system has been dis- 
turbed so that the instantaneous temperature, pressure, etc. can be defined via flS]) for the 
system even if they are changing with time. In other words, there is partial equilibrium in 
the isolated system. We then apply this thermodynamics to study glasses and clarify the 
concept of the Active temperature T n (t) widely used in the study of glasses by identifying it 
as a thermodynamic quantity; see fl25|) . Our analysis shows that the Active temperature has 
the same value even if we change the relaxing quantity. This temperature is not identical 
to but is related to the instantaneous temperature T(t) in a glass; see (1261) . We use this 
relationship to establish the Tool-Narayanaswamy equation (12 7p for the relaxation time on 
a solid theoretical ground. This form does not change even if we have more than one kind of 
slow relaxation; the latter results in many different Active temperatures, one for each kind 
of slow relaxation. 

We should finally contrast our approach with other approaches available in the literature. 
We can use (12 2 aft to express dS(t) = dS c (t) + dS n (t). However, each term is still multiplied by 
T(t), implying that dS e (t) or dS n (t) are not multiplied by their respective temperatures T or 
T n (t). Thus, one cannot consider S e (t) and S n (t) as separate in the first law ffTO]) or in fT2~Tl) . 
This should be contrasted with the approach developed in 22J, where the entropy is divided 



into fast and slow dof; see also [23j]. More recently, Moller et al 



241 ] have used the approach 



of de Donder to study glasses by considering a single structural order parameter; however, 
the concept of the Active temperature was not analyzed. A recent approach by Wolynes 
(b)] provides a local description of the relaxation in an inhomogeneous mosaic form, but the 
interest is in the dynamics, whereas our focus is not on any particular dynamics. 
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